International Journal of Power Electronics and Drive Systems (IJPEDS) 
Vol. 13, No. 3, September 2022, pp. 1486~1497 
ISSN: 2088-8694, DOI: 10.1159 1/ijpeds.v 13.13.pp1486-1497 O 1486 


Minimization of total harmonic distortion in neutral point 
clamped multilevel inverter using grey wolf optimizer 


Fahmi Ahyar Izzaqi, Novie Ayub Windarko, Ony Asrarul Qudsi 


Electrical Engineering Department, Politeknik Elektronika Negeri Surabaya, Surabaya, Indonesia 


Article Info 


ABSTRACT 


Article history: 


Received Mar 18, 2022 
Revised May 22, 2022 
Accepted Jun 15, 2022 


The inverter has been attracting researchers for their application in 
renewable energy. So far, multilevel inverter is considered as low distortion 
class, which produces multilevel output voltage imitating a pure sine 
waveform. However, the needs for free distortion of output voltage have 
been motivating to improve multilevel pulse width modulation PWM 


generation method. In this paper, the modified PWM technique is proposed 
to reduce the voltage total harmonics distortion (THD) of multilevel inverter. 
Keywords: This modulation technique is then applied to control a single-phase three- 
level neutral point clamped multilevel inverter (NPC-MLI. Grey wolf 
PA optimizer (GWO) algorithm is utilized to generate optimal amplitude and 
Modified P WM phase shift of modified reference signal. The GWO algorithm is then 
Multilevel inverter compared with other optimization algorithms such as differential evolution 
Optimization algorithm (DE), human psychology optimization (HPO), and particle swarm 
Total harmonic distortion optimization (PSO) to evaluate their performance in harmonic minimization. 
The performance of the proposed work is validated through simulation and 
experimentation on a prototype. The results show that the modified PWM 
technique optimized with GWO can reduce THD on NPC-MLI output 
voltage. 


Grey wolf optimizer 


This is an open access article under the CC BY-SA license. 


Corresponding Author: 


Fahmi Ahyar Izzaqi 

Electrical Engineering Department, Politeknik Elektronika Negeri Surabaya 
ITS Raya st, Keputih, Sukolilo, Surabaya, East Java 60111, Indonesia 
Email: fahmizzaqi@pe.student.pens.ac.id 


1. INTRODUCTION 

Since several decades ago, multilevel inverter multilevel inverter (MLI) has been attracting to 
convert direct current DC voltage to alternating current AC voltage in renewable energy applications [1]. The 
MLI are widely used as voltage converters from renewable energy source such as solar panels, wind turbines 
and others, to supply electricity to AC loads in stand-alone systems [2] or to the grid-connected [3]. High 
efficiency power conversion systems are expected to deliver high performance by reducing conduction losses 
in the switching process and reducing total harmonics distortion total harmonics distortion (THD) [4]. So far, 
MLI is well known has low voltage distortion and can operate at high power and high switching frequency 
[5]. Under ideal condition, the MLI should provide pure sine wave for electric loads. The MLI can produce 
staircase output voltage which close to a sinewave form. However, the MLI’s staircase output voltage is non 
ideal sine wave and still contains harmonics due to the switching pattern generation [6]. To overcome this 
issue, many type of MLI has been developed includes cascaded H-Bridge MLI, flying capacitor MLI, and 
neutral point clamped multilevel inverter neutral point clamped multilevel inverter (NPC-MLI) [7]. Due to 
MLI circuit complexity, several researchers have proposed to reduce the number of MLI’s components [8]. 
Nevertheless, the NPC MLI has advantages of easy control and low cost hardware implementation among the 
MLIs [9]. 
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Beside the topological solutions, reducing THD of output voltage by generating proper switching 
pattern modulation is also a hot issue of research in MLIs [10], such as multi-carrier PWM [11], selective 
harmonic elimination pulse width modulation PWM [12], and space vector modulation [13]. Many studies 
have been performed to reduce the harmonics of the output voltage of a MLI [14]. Several modulation 
techniques have been reported to selectively eliminate lower order harmonics of MLIs staircase output 
voltage [15]. Moreover, injection techniques are also common by injecting various reference signals to 
modulated signal [16]. Many researchers have proposed new injection techniques such as third harmonic 
injection PWM [17], third harmonic injection space vector PWM [18], and generalized discontinuous PWM 
[19]. Each injection technique has shown specific and various THD value of the output voltage. Recently, 
several modulation techniques to minimize the THD of MLI output voltage was improved by utilizing 
metaheuristic algorithm [20]. Improved whale optimization algorithm is used to optimize the switching 
angles of a single phase 5 and 7 level multilevel inverters (MLIs) [21]. Grasshopper optimization algorithm is 
used for optimizing switching angles applied to a three-phase 7-level cascaded H-bridge multilevel inverter 
[22]. To minimize the THD in cascaded MLI, the hyper-spherical search algorithm is used [23]. The use of 
metaheuristic algorithm to optimize injection techniques in NPC-MLI is proposed, such as genetic algorithm 
(GA) [24] and differential evolution (DE) [25]. However, drawbacks in metaheuristic algorithms can restrict 
the minimization of THD. For example, insufficient local search ability, slow convergence rate and 
premature convergence are main limitations associated with GA [26]. Also, the DE performance are fully 
influenced by the parameter of scaling factor and the crossover rate [27]. So, the determination of those 
control parameters is significant in the mutation and crossover operations. Hence, to improve the THD 
minimization, researchers have been promoting novel modulation and higher performance optimizer algorithms. 

Based on the above mentioned, this paper is intensively proposed a low-cost approach to reduce 
THD by injection technique to modulated signal for NPC-MLI. This paper also provides discussion and 
analysis of the injected PWM optimization using grey wolf optimizer GWO. GWO is a meta-heuristic 
method that is guaranteed to be very efficient and effective when applied to solve constrained problem or 
unconstrained problem, fast to achieve convergence, and avoid local optimization [28]. GWO imitates the 
hunting strategy and social hierarchy of grey wolves. GWO has four hunting strategies of searching, 
encircling, hunting, and attacking to achieve the prey location. The social hierarchy consists of alpha, beta, 
delta, and omega wolves. Alpha, beta, and delta wolve means the first, second and third best solution for 
optimization respectively. Omega wolves means worse solutions than delta wolve. GWO is used to find the 
best injected signal parameters of frequency, amplitude, and phase shift. The PWM modulation techniques is 
applied to NPC-MLI, then tested in PSIM software simulation and verified under experimental works. 


2. RESEARCH METHOD 
2.1. Neutral point clamped—multilevel inverter 

The inverter is a device that converts a DC voltage into AC voltage. Conventional inverters 
typically have three output voltage levels, +Vdc, -Vdc, and zero. Multilevel inverter is a type of inverter that 
can generate multilevel voltage and current levels. NPC-MLI also known as diode clamped multilevel 
inverter (DC-MLI). The NPC-MLI consists of two pairs of switches in series (upper and lower), parallel to 
two capacitors in series where the anode of the upper diode is connected to the midpoint of the capacitor and 
the cathode to the midpoint of the top half diodes of the switch pair. The cathode of the lower diode is 
connected to the mid-point of the capacitor which enables the main DC voltage to be divided into smaller 
voltages (1/2 VDC). The mid-point of the two capacitors can be defined as the “neutral point” [24]. Figure 1 
shown a schematic of single phase three level neutral point clamped multilevel inverter. The NPC-MLI has 
(m-1) the number of DC capacitors, (m-1)*(m-2) the number of clamping diodes, and 2*(m-1) the number of 
switching components per phase, m is the level of the multilevel inverter. 


2.2. Modified PWM technique 

The switching process is the cause of the THD of the NPC-MLI output voltage. As a result, an 
appropriate modulation technique is required to reduce the effect of harmonic distortion. In this paper, a 
modified PWM technique is proposed. The conventional modulation technique that is commonly used is 
sinusoidal pulse width modulation SPWM. This modulation technique is formed by comparing the sine 
reference signal with the triangular carrier signal. NPC-MLI requires a few semiconductor components so 
that more than one reference signal and a carrier signal are required. 

The modified PWM techniques can be observed in Figure 2. The modified reference signal is 
produced by combining two signals of reference signal and injected signal. The reference signal is sinus 
waveform with fundamental frequency, and the injected signal has designated value of frequency, phase shift, and 
amplitude. The injected signal can be single or multiple signals. The mathematical equations for signals are as: 
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V,(t) = A, Sin2rf,t (1) 
V; (t) = per An (Sin2nf,t + On) (2) 


Which V, is reference signal and V; are injected signals. 
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Figure 1. Schematic of single phase three level NPC-MLI 
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Figure 2. Schematic of the modified 


In (1) is a sinusoidal reference one that has an amplitude of 1 and fundamental frequency of 50 Hz. 
in (2) is signal of injected signal. The injected signal aims to modify the reference waveform. This signal has 
the value of amplitude (A,,), frequency (fn) and phase shift (6,,), which their best value will be found through 
the Grey Wolf Optimizer tracking process. The two signals are added, resulting in the mathematical in (3), 
and for example in (4). 


Ving(t) = Ay Sin2mf,t + Zz- An (Sin2mfyt + On) (3) 
Ving(t) = Sin2n50t + 0.1 {Sin2n750t + 180°} (4) 


Which Vms is modified sinusoidal reference signal. 

Figure 3 illustrates the process of generating a modified PWM signal using in 1-3. Figure 3(a) shows 
a sine reference wave one based on in (1), with an amplitude of 1 and a fundamental frequency of 50 Hz. 
Also, Figure 3(b) shows an injected signal waveform based on in (2) which has an amplitude of 0.1, a 
frequency of 750 Hz, and a phase shift of 180 degrees. Figure 3(c) shows both signals are combined and 
generate a modified reference signal that will be compared to the triangular carrier signal and form a 
modified PWM. 


2.3. Grey wolf optimizer for minimization THD 

The grey wolf optimizer (GWO) is a meta-heuristic method for finding optimal solutions of 
numerical problems. This method was adapted from the hunting behaviour of a herd of grey wolves. A wolf 
pack has a division of the leadership hierarchy, namely alpha (a), beta (B), delta (6) and omega (œ). Alpha is 
the best solution and as the leader of the herd. Beta and delta are the second and the third best solution that 
stands for alpha. While omega is the rest of the herd [20]. Figure 4(a) shows the hunting behaviour of grey 
wolves and Figure 4(b) shows flowchart of the GWO algorithm. 
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Figure 3. Modified reference signal generation scheme: (a) sinusoidal reference, (b) injected signal, 
and (c) modified reference signal 
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Figure 4. Illustration of GWO algorithm: (a) the hunting behavior of grey wolves [28] 
and (b) GWO flowchart 
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Based on the behaviour of wolf hunting, this method consists of 4 steps for obtaining the best 
solution. These steps include tracking, encircling, hunting, and attacking. It has been adapted as a method to 
find the best parameters in the modified PWM modulation process. The process of this algorithm can be 
explained as: 

— Searching 

Searching is the first step for a pack of wolves to hunt prey. In this step all members of the herd will 
spread out in a certain area to find prey. In its application to modified PWM, this step is carried out by giving 
the initial value of the solution or the initial position of the herd randomly with a certain range of the parameters of 
frequency, amplitude, and phase shift. After spreads, each member's position in this herd will be evaluated. 

— Encircling 

Encircling can be illustrated as a step to surround prey. Each member of the flock will keep their 
distance and move progressively closer to the prey. This siege will be led by alpha (a) and assisted by beta (B) and 
delta (5). While the movement of omega will follow alpha, beta and delta. This step is adapted into the following. 


= 


=|¢-X,(t) - p | (5) 


SI 


X(t+1) =X,(t)-A-D (6) 


In the equation, t is the current iteration value. While D, A, and C are vector coefficients. X, is the position 
vector of the prey and X is the position vector of the wolf. Vector values A and C can be calculated using the 
following equation. 


Ā=2ä -å (7) 
Č=2:7 (8) 


The value of a decreases linearly from 2 to 0 according to the iteration value. While the values of 7, and rz 
are generated randomly. 
— Hunting 

In an abstract search space, the optimal location (prey) is unknown. Alpha is the best candidate 
solution, whereas beta and delta have better insight about potential prey locations. The first three best 
solutions obtained should update their positions accordingly to the best search agent positions. The following 
formula is proposed. 


Da = |G.Xq — X|,Dg = |02-Xg — X|, Ds = |c3.Xs — X (9) 

X, = Xa — Aı.Da» Xz = Xg — Az. Dg, X3 = X6 = Az. Dy (10) 

> i+ +ž. 

X(t+1) = ae (11) 
— Attacking 


Attacking the prey is the final step in the wolf hunting process. The attacking process is carried out 
when the prey has stopped moving. In its application to modified PWM optimization, this process is carried 
out when the solution results from each member of the herd has reached the convergence value. This value 
will then be performed as a best solution. Illustration of the grey wolf hunting process shown and the 
flowchart of the GWO algorithm is shown in Figure 4. 


3. SIMULATION AND EXPERIMENTAL RESULTS 
3.1. Simulation results 

PSIM software with simplified C-Block are used to simulate the modified PWM technique 
optimized with the algorithms. According to Figure 1, the parameters of the NPC-MLI used in the simulation 
are shown in Table 1. To simplify the optimization, the injected signal is single waveform, and the frequency 
is determined manually. The frequency values are determined up to 21* harmonics. They are chosen to 
explore the best result from various value of injected frequency. Then, the amplitude and phase shift of the 
injected signal is optimized by GWO. Figure 5 shows an example of the tracking process results of injected 
signal with frequency of 650Hz. It can be seen from Figures 5(a) and 5(b), amplitude and phase shift injected 
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are searched initiated from zero then approached to best value. Figures 5(c) and 5(d) show the THD results 


from injected amplitude and phase shift. At time of around 6.1s, the best value of amplitude and phase shift 
are convergence. 


Table 1. NPC-MLI simulation parameters 


Components Value 
DC Source 311 V 
Capacitor 1 270uF 
Capacitor 2 Load 270uF 100 Ohm 
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Figure 5. Simulation result in frequency 650Hz of optimized injected signal, (a) amplitude, (b) phase shift, 
(c) THD, and (d) the best THD reached 


Optimization is performed to 10 various frequencies of injected signals. GWO is tested and 
compared its performance to other optimization methods, such as differential evolution (DE), human 
psychology optimization (HPO), and particle swarm optimization (PSO). Figure 6 shows comparison charts 
of the optimization results of the GWO, DE, HPO, and PSO algorithms. The charts show the best amplitude, 
the best phase shift, convergence time and the best THD. Convergence time in those charts mean the time 
needed to reach the best THD, such as shown in Figure 5(d). Convergence time is shown to represent the 
speed of optimization process. Figures 6(a) and 6(b) show the results of amplitude and phase shift. 
Figure 6(b) show clearly that GWO can find the best value of phase shift by avoiding trapped in local 
maxima. Contrarily, other algorithms trapped in local optima which indicated by extreme phase shift value of 
0° or 180°. Consider to Figure 6(c), the lowest average THD is reached by GWO at 28.5%, then HPO 28.6%, 
DE 28.7% and PSO 28.8%. Figure 6(d) shows that DE is the slowest to convergence. Furthermore, the fastest 
average convergence time is reached by GWO at 6.079s then PSO 6.080s, HPO 6.164s and DE 37.118s. 
Figure 7 shows comparison of the THD minimization results plotted by time using GWO, DE, HPO and PSO 
at a frequency of 650Hz of injected signal. The injected signal resulted from GWO yielding a THD of 
26.44%. The GWO process runs and finds the best solution in 6.14s. HPO yield THD of 26.61%, DE 26.62% 
and PSO 26.96%. HPO yield convergence time of 6.18s, PSO 6.44s and DE 37.8s. 

Figures 8(a) and 8(b) shows the comparison of NPC-MLI output voltage waveform using SPWM 
and modified PWM. Since both of output voltage has same amplitude modulation, the injected signal effects 
to modified PWM can be identified by comparing the different pattern Figure 8(a) and 8(b). Figure 8(c) show 
clearly the THD level of modified PWM is lower than SPWM. Table 2 provides a comparison of the THD 
level of the NPC-MLI output voltage in percent. The SPWM has a THD of 31.42% and the modified PWM 
has a THD of 26.44%. Figure 9(a) experimental setup of proposed work and 9(b) hardware of NPC-MLI. 
Figure 10(a) harmonics spectrum of SPWM and 10(b) harmonics spectrum of modified PWM. 
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Figure 6. Optimization results from various algorithms, (a) amplitude (b) phase shift (c) THD, 
and (d) convergence time 
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Figure 7. Comparison of the THD minimization results using GWO, DE, HPO and PSO at a frequency of 
650 Hz of injected signal 
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Figure 8. Output voltage of NPC-MLI, (a) SPWM, (b) modified PWM, and (c) SPWM THD and modified 
PWM THD on NPC-MLI in simulation 


Table 2. THD comparison of modulation techniques in percent using simulation 
Modulation Technique _ THD (%) 
SPWM 31.42 
Modified PWM 26.44 
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Figure 9. The result of the proposed work in hardware: (a) experimental setup and (b) NPC-MLI prototype 


Figure 10. Harmonics spectrum: (a) SPWM and (b) modified PWM 


3.2. Experimental results 


To verify the proposed modulation technique and to validate the simulation results achieved by 
GWO algorithm, a single-phase NPC-MLI prototype was built. Table 3 shows the components used to 
develop NPC-MLI. Figure 9(a) shows the setup of experimental system and Figure 9(b) shows the zoomed 
picture of NPC-MLI. The Nucleo STM32F446RE microcontroller is used to generate pulses for NPC-MLI 
switching. Picoscope is used to measure and to analyze the waveform, harmonics spectrum and THD. 
Programmable DC Power supply TDK-Lambda is used to provide constant DC voltage input of NPC-MLI. 
The NPC MLI is tested with an DC input voltage of 311V and is loaded with a resistive load. 


Table 3. Components used for experiment 


Components Specifications 
IGBT Module K40H1203 (Infineon) 
Microcontroller Nucleo STM32F446RE 
Optocoupler FOD3182 (Fairchild) 
Fast Diode MUR 1560 
Gate Driver Power Supply DC 12V, HLK-PM12 
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NPC-MLI is operated to verify the modulation of SPWM and modified PWM technique with the 
best parameters from the simulation results. Both of SPWM and the modified PWM is set to amplitude 
modulation of 1. Furthermore, the modified PWM is injected by signal of 650Hz with an amplitude of 
0.07519 and phase shift of 180 degrees. Figures 11(a) and 11(b) show the comparison results of experiments 
between output voltage of SPWM and modified PWM. The results are similar to simulation results as shown 
in Figure 8. The output voltage is recorded to determine the THD. Figure 12 shows the spectrum and THD 
results. Both of modulation produce harmonic spectrum at very low level and at high frequency. The high 
harmonics of frequency spectrum can be advantage for the suppression of distortion by only using small LC 
filter. From Figure 12, SPWM produces higher THD than modified PWM. Figure 12(a) shows the THD 
measurement using the spectrum mode of picoscope 4824 is 28.15%. Figure 12(b) shows the output 
waveform of a NPC-MLI using Modified PWM switching. The THD is noticeably less than SPWM. 
Figure 12(b) The measured THD is 26.4% which proves that switching Modified PWM optimized by GWO 
produced lower THD than SPWM. 
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Figure 11. Output voltage of NPC-MLI: (a) SPWM and (b) modified PWM 
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Figure 12. THD measurement using spectrum mode of picoscope 4824: (a) THD of SPWM and 
(b) THD of modified PWM on prototype system 


4. CONCLUSION 

This paper presents modified PWM with optimization using the GWO algorithm as an alternative 
modulation technique to reduce voltage harmonics in NPC-MLI. The values of frequency, amplitude, and 
shift angle are optimized to produce low THD. The GWO algorithm is then compared to other algorithms 
such as DE, HPO and PSO to test their performance in harmonic reduction and convergence level. As a 
result, GWO can optimize modified PWM to obtain global optima with a minimum. To validate the proposed 
theory, a prototype has been built and tested to compare between common SPWM and Modified PWM 
optimized by GWO. The lower THD produced by proposed method can be applied to increase the 
performance of NPC-MLI. Furthermore, the low THD NPC-MLI can help to promote the next level of 
renewable energy applications. 
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